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ABSTRACT 


The  multiple  target  ocean  surveillance  contact  correlation  problem  can  be  decoupled  roughly  into  two  parts.  The 
first,  extensively  treated  in  the  literature,  involves  gcolocation  information  only,  and  is  normally  analyzed  through 
use  of  a  bank  of  Kalman  filters.  The  second  is  concerned  with  non-geolocac lor  attribute  information.  Typically,  the 
letter  includes  all  data  obtained  through  linguistic,  visual,  or  discrete  valued  numerical  sources. 

In  this  paper,  a  procedure  is  proposed  which  analyzes  non-geolocat ton  attribute  information  and  yicls  a  posterior 
possibility  distribution  of  target  correlations.  In  torn,  this  result  raa>  be  utilized  to  compute  the  overall  posterior 
probability  distribution  of  correlations.  The  key  factor  in  this  procedure  is  a  theorem  which  shows  that  a  uniform/ 
most  accurate  confidence  sot  exists  which  is  determined  by  a  single  possibility  distribution  chat  is  feasible  to 
compute. 


INTRODUCTION 

The  theme  of  the  48ch  MORS,  "Military  Operations 
Research  Techniques  for  the  80's",  is  bcjjig  addressed  }.i\ 
many  different  ways  throughout  the  symposium.  Most  of 
the  approaches  are  either  deterministic  or  probabilistic 
in  nature,  reflecting  the  trend  in  current  analysis,  yet 
there  are  many  problems,  military  and  non-military,  that 
can  not  be  easily  formulated  in  probabilistic  (nor  deter¬ 
ministic)  terms.  For  example,  when  information  is  gleaned 
from  strictly  human  operator  sources  such  as  chrough  vis¬ 
ual  sitings  or  judgments  based  on  experience,  linguistic 
discretions  may  be  chc  prime  data  base.  Such  information, 
although  often  ambiguous  or  vague,  certainly  may  prove 
useful  -  often  in  conjuncion  with  "harder"  statistical 
information. 

In  the  39ch  MORS,  Dr.  J.T.  Dockery  most  aptly  showed 
how  fuzzy  sec  techniques  could  be  used  in  military  prob¬ 
lems  (1).  Among  ocher  papers  demonstrating  the  feasibility 
of  fuzzy  sec  theory  and  techniques  in  a  military  context, 
one  may  include  chc  basic  work  of  Watson  cc  al.  (2)  and 
this  author  ((3),  Examples  2  and  4).  This  paper,  coo, 
will  deal  with  one  aspect  of  a  class  rf  military  problems 
-  the  multiple  target,  multiple  sensor  contact  correla¬ 
tion  problems  -  chrough  the  use  of  fuzzy  set  theory. 

The  tneory  of  fuzzy  secs  is  relatively  new,  being 
formulated  fully  for  the  first  time  by  Zadeh  in  1965  (4). 
(For  earlier  attempts,  see  for  example,  the  interesting 
work  (1937)  of  Shirai  (5)  and  (independently)  Black  (6); 
Sheppard's  psychological  quantifications  (7);  Watanabo’s 
contributions  (aimed  coward  quantum  mechanics) (8) ;  and 
Klaua's  work  based  on  many-valued  set  theory  (9).)  At 
present,  much  work  has  been  done  justifying  rigorously 
chc  Important  role  that  fuzzy  sec  theory  plays  with 


respect  to  probability  cheory((lO),  (11)).  (Future  work 
will  address  the  fuzzy  set  approach  to  data  sampling 
and  the  Laws  of  Large  Numbers  and  Central  Limit  Theorems) 
Essentially,  it  can  be  shown-emphasiziug  the  connection 
of  fuzzy  sec  modeling  to  vagueness  of  information-  that 
any  fuzzy  set  can  be  identified  in  a  natural  way  with 
an  equivalence  class  of  random  secs  (in  general,  many, 
unless  the  fuzzy  sec  is  an  ordinary  set)  and  chat  fuzzy 
and  random  set  operations  correspond.  Nevertheless, 
many  of  the  difficulties  inherent  in  chc  probabilistic 
approach  can  be  avoided  by  use  of  the  well-developed 
calculus  of  fuzzy  set  operations  and  relations.  Specific¬ 
ally,  determination  of  joint  distributions  and  integra¬ 
tion  of  functions,  a  necessary  factor  iai^robabHiscic 
modeling,  is  replaced  by  much  simpler  operations.  (Sec, 
for  example,  Dubois  and  Prade's  comprehensive  text  (12), 
illustrating  these  procedures.) 

On  chc  other  hand,  it  is  not  the  intent  of  chc  above 
exposition  to  claim  that  fuzzy  set  techniques  should  re¬ 
place  probabilistic  ones.  Rather,  it  is  chat  Chc  two 
approaches  to  modeling  and  manipulating  uncertainties 
may  be  used  compatibly  (sec  (13)):  when  numerical  des¬ 
criptions  are  available  with  well  defined  probability 
distributions  present,  use  a  probabilistic  approach; 
when  linguistic  descriptions  or  ocher  vague  information 
is  present,  use  a  fuzzy  sec  -  or,  to  establish  an  anala- 
gous  terminology  -  posslbiliscic  approach. 

Although  an  elementary  treatment  of  possibility 
theory  and  techniques  will  not  be  given  here  (see,  for 
example,  (12)  or  (l^)).  one  basic  property  contrasting 
with  probability  theory  must  be  presented.  In  probabil¬ 
ity  theory,  the  well-known  concept  of  a  probability  dis¬ 
tribution  plays  a  key  role.  In  possibility  chcory,  the 
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possibility  distribution  plays  a  key  role.  In  the  latter, 
equivalently  described  by  a  fuzzy  set  membership  function 
generalizing  the  concept  of  the  membership  or  character¬ 
istic  function  of  an  ordinary  set  (by  allowing  values 
that  nay  often  be  neither  zero  nor  one,  but  intermediate), 
the  sun  (or  integral,  if  appropriate)  of  the  possible 
values  need  not  add  up  to  one  -  as  opposed  to  the  sun  of 
values  described  by  a  probability  distribution.  This 
concept  will  be  clarified  further  through  the  procedure 
developed  later  in  the  paper. 

Consider  now  the  general  contact  correlation  problem. 
This  nilitary  problen  has  had  a  long  history,  beginning 
with  the  early  fornulations  of  Sittler  (196G) (15)  and 
Wax  (1955) (16)  down  through  the  present.  In  1979,  well 
over  300  papers  were  gathered  at  the  Naval  Research  Lab¬ 
oratory  as  part  of  the  Naval  Ocean-Surveillance  Corre¬ 
lation  Handbook  Project.  (Sec  (17)  for  a  conpendiun, 
overview,  and  analysis  of  raany  of  these  papers  in  the 
correlation  field.  See  also  (18)  and  (19)  for  further 
analysis  of  the  general  problen.)  Basically,  the  problen 
can  be  stated  as  a  data  partitioning  one,  where  that 
partitioning  iu  sought  where  each  conponent  represents 
sensor  or  huc^n  source-gathered  information  (over  some 
sampling  tine  period)  partaining  to  the  sane  target  or 
object.  The  information  nay  be  roughly  divided  into  three 
classes:  1,  geolocation  sensor-obtained,  such  as  bearing 
and  range  measurements;  2,  false  alarm,  an  encompassing 
term  used  to  denote  data  arising  from  sources  which  are 
either  of  no  interest-icebergs,  neutral  ships  -  or  false 
signals  due  to  reflection  and  scattering,  for  exanplc; 
and  3,  attribute  information.  Models  which  address  the 
first  two  types  of  information  dominate  the  correlation 
field.  (Again,  see  (17)  for  further  details.)  However, 
the  work  of  Reid  (20)  and  Bowman  (21),  among  others,  in 
attempting  to  incorporate  the  third  category  of  data 
should  be  cited.  Nevertheless,  to  the  author's  knowledge, 
no  approach  to  chirproblen  through  possibility  clrcoryr- 
until  now,  has  been  undertaken.  (Figure  1  illustrates  a 
typical  correlation  situation.) 

UNCLASSIFIED 
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(UlTypical  Correlation  Problem 
Figure  1 


What  is  meant  by  attributes  in  the  context  of  the 
contact  correlation  problem?  An  example  should  suffice: 
Target  A  may  have  associated  with  it  maneuvering 
characteristics,  signal  frequency  information,  and  a 
visual  siting  indicating  an  irregular  design.  Target  B 
may  also  have  similar  maneuvering  characteristics,  signal 
pattern,  and  a  tentative  classification.  In  addition, 
both  targets  have  been  related  by  intelligence  Informa¬ 
tion,  although  certain  discrepancies  appear.  Finally, 
Target  A  has  been  determined  to  have  on-board  three  radar 
sensors,  while  Target  B  is  known  to  have  "several  sensor 
systems  operating,  and  appears  to  be  heading  to  port". 
Should  Targets  A  and  B  be  correlated,  tentatively  cor¬ 
related  until  statistical  data  (gcolocacional)  is  avail¬ 
able,  or  should  the  two  be  ruled  out  for  possible  corre¬ 
lation?  Another  example  could  be  generated,  where  a  mix 
of  statistical  and  verbal  descriptions  are  present  for 
both  targets  in  question.  In  any  case,  certainly, 
wherever  possible,  the  experience  of  human  decision 
makers  should  also  be  taken  into  account  -  as  well  as 
available  automatic  correlation  algorithms  -  in  order 
to  utilize  fully  the  information.  (See  figure  2.) 
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(U) Example  of  Attributes 
Figure  2 


Following  simplifying  assumptions  concerning  the 
statistical  dependencies  of  the  relevant  variables  and 
the  accuracy  of  the  overall  modeling  (which  may  change 
drastically  over  a  sufficient  period  of  time),  the 
following  theorem  concisely  shows  the  role  that  attribute 
probabilities  play: 

Theorem  Structural  decomposition  of  correlation  problem. 

Let  denote  any  partitioning  of  data 

Z*(2g,Zafz^>  up  to  time  tj  ,  where  the  subscript  f  indi¬ 
cates  false  alarm,  g  Indicates  gcolocacional,  a  indicates 
attributes,  and  +  denotes  the  collection  of  all  compo¬ 
nents  of  Q  which  correspond  to  targets  of  interest.  Then 
assuming  statistical  dependencies  only  occurring  in  the 
conditional  forms  remaining, 

P(Q'Z)  -  p(ZgUa,Q).p(Z£|Qf).p«)l7.;>)/P(Zg.ZflZ3)-  <*) 
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(The  proof  of  the  theorem  follows  directly  froa  che  prop¬ 
erties  of  conditional  probabilities.) 

■ 

See  (17)  and  (18)  for  background  and  related  results. 
Note  that  the  theorea  states  that  the  posterior  probab¬ 
ility  distribution  of  the  data  partitionings  depends 
directly  on  the  posterior  distribution  of  che  data  par¬ 
titionings  given  the  attributes.  The  first  factor  in 
eq.  (1)  corresponds  to  the  Hainan  filter  innovations 
(under  Che  standard  assuaptions  of  Causs-Markov  linear 
target  state  relations  and  measurements),  while  che 
second  factor  nay  be  obtained  in  several  different  ways, 
according  to  the  nodel  assumed  for  Che  dispersion  and 
occurence  of  false  alams.  The  divisor  is  a  function  of 
Z  only  and  not  of  Q.  Hence,  this  tern  plays  no  role  when 
optinun  Q  is  sought,  i.c.  that  value  of  Q  -  subject  to 
feasible  search  constraints  -  which  maximizes  the  pos¬ 
terior  probability  p(QlZ).  Finally,  it  should  be  noted 
that  the  above  theorem  is  concerned  solely  with  probab¬ 
ilities,  not  possibilities.  However,  che  renainder  of 
this  paper  will  be  devoted  to  showing  how  possibilities 
(and  fuzzy  set  theory)  can  be  used  to  obtain  the  desired 
probabilities. 
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(UjExaaplc  1  of  a  Posterior  Possibility 
Distribution  for  an  Attribute 
Figure  3 

III.  Consider  nou  tae  set  of  experts  y-.y,,...  Eacli 
expert  y^  Is  asked  the  following  question: 


The  naln  goal  of  chls  paper  can  nou  be  stated:  To 
shoo  chat  there  exists  a  feasible  and  cacheaatlcally 
justifiable  procedure  for  obtaining  the  posterior  possi¬ 
bility  distribution  of  q  given  Za  ,  and  in  turn  using 
this  expression  to  generate  a  naturally  corresponding 
evaluation  for  p(Q|Za),  which  froa  eq.  (1)  aay  be  used 
in  determining  the  overall  value  of  p(QiZ). 

BASIC  MODEL 

The  procedure  for  obtaining  che  posterior  possibility 
distribution  of  data  partitionings  given  the  relevant 
attribute  information,  can  be  conveniently  divided  into 
eight  steps. 

’  -1  -  ’  —  —  *  >  /  .  . 

I.  A  taxonomy  of  accrlbu  c  relations  is  developed  us¬ 
ing  heuristic  or  statistical  procedures.  (Sec  for 
example,  che  approach  of  Novakovska  (22).  The  goal  here 
is  che  definition  of  a  sec  of  relatively  independent 
critical  attributes  for  the  correlation  problca  at  hand. 
For  example,  this  list  could  include  A|-  degree  of 
observed  maneuvering,  Aj-  signal  frequency  character¬ 
istics  (this  could  be  further  subdived),  A}«  bearing 
information,  classification  (although  this  attribute 
nay  depend  to  a  lafge  degree  on  core  primary  ones) , 

A^«  nusber  of  sensor  systens  on-board,  target  ident¬ 
ification  nuaber.  Ay-  visual  characteristics  (again,  sub- 
ivision  into  core  specific  categories  is  norc  ceaning- 
ful).  Call  the  final  set  of  primary  attributes 

C  Aj,A2,..  .Aj^ 

II.  The  domain  of  possible  values  or  confusables  that 
each  attribute  can  assume  is  determined.  In  order  to 
accomplish  this  and  related  tasks,  a  panel  of  experts 
nust  be  available  for  querying.  The  objective  in  this 
phase  is  che  nodeling  of  the  posterior  attribute  possi¬ 
bility  distributions  given  typically  by 


^(Aj|zj)  :  Dj  —  I0-1!  ' 


j  -  1,..,M,  where  each  Aj  is  an  attribute,  Z  j  <■  Dj  is 

arbitrary  representing  any  potential  observed  (data) 
value  of  Aj  ,  Dj  being  the  domain  or  set  of  all  possible 
values  of  Aj,  usually  seasured  in  some  convenient  dimen¬ 
sional  units,  (Sec  figure  3.) 


"What  is  the  possibility  that  given  you  have  observed 
data  Zj  as  a  value  for  attribute  Aj,  that,  say,  Vj  is  the 
actual  (or  equivalently,  a  confusabie)  vilue  for 
attribute  Aj?" 


Each  value  Vj  is  then  coved  around  freely  in  domain 
Dj,  and  in  turn,  the  Zj  *s  arc  moved  about  in  Dj,  for 
each  j,  j-l ,2, . . ,M.  As  a  check,  when  Vj-Zj,  the  response 
possibility  should  be  unity  or  nearly  so.  Thus,  symbol¬ 
ically  <p  _  (V.)  represents  the  possibility 


that  Vj  has  attribute  Aj.  given  Zj  is  observed  by  expert 

yfc.  3ayes'  cheorca  {see  phase  VI)  in  a  fuzzy  context 
could  have  been  used  here,  but  would  entail  prior  distri¬ 
butional  assumptions  as  well  as  additional  calculations. 
Instead,  the  posterior  possibility  distributions  arc 
obtained  divectly.  (Sec  also  figure  4.) 

\ 
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Figure  4 
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IV.  a  simple  calculation  shows  chat  if  s<  is  the 
number  oi  elements  in  domain  D j ,  then  (Sj)2  calculations 
are  required  in  step  III  for  each  corresponding  poster¬ 
ior  possibility  distribution.  One  way  of  avoiding  ail  of 
these  calculations  is  to  use  *  where  appropriate,  of 
course,  such  as  for  attributes  Aj  and  A3  in  the  example 
in  step  I  -  a  translation  parameter  family  of  functions. 
Thus,  some  of  Che  possibility  distributions  obtained  in 
III  may  be  put  in  the  form 

^(Ajl  Zj.yk)  <V  -  P<Aj|  yk)  <V2J>.  (3) 

V.  The  next  step  in  the  procedure  is  the  determination 
of  the  sec  of  relevant  rules  used  in  the  decision  pro¬ 
cedure  by  the  experts  when  correlation  is  to  be  carried 
out  based  on  the  available  information.  More  specifi¬ 
cally,  the  rules  considered  here  arc  of  the  form 


Other  rules  cay  cake  Into  account  geographical 
barriers  or  physical  constraints. 

Care  must  be  taken  here  in  the  modeling  of  each  im¬ 
plication  (or  modus  ponensk  rule  Ej  chat  the  implication 
is  based  on  an  ideal  situation,  i.c. ,  no  error  is  assumed 
for  the  observation  of  the  attribute  values;  errors  are 
accounted  for  separately  (and  then  combined  optimally 
with  the  rules  in  step  VIII).  Although  the  model  can 
handle  match-no  match  situations  in  a  non-trivial  manner, 
analagous  to  the  testing  of  hypotheses  in  a  probabil¬ 
istic-statistical  situation,  more  flexibility  is 
achieved  when  the  catch-no  match  situations  are  replaced 
by  fuzzy  gates  (analagous  to  the  introduction  of  random¬ 
ness  for  the  parameters  involved  in  the  testing  of  hypo¬ 
theses  analogy). 


If  attributes  As  and  As.  and.,  or  A,  or  A.  or.. 


and 


*2 


or 


etc.  are  present  at  accual  values  Z 


Ji- 


Zj^,  Zjj,...,  respectively,  then  correlations 

or  equivalently  data  partitionings  Q|,  • _  are  poss¬ 

ible  utch  Qj  being  cost  likely,  Q2  less  likely,  etc." 


In  tens  of  attributes  or  fatty  sets  the  above  ex¬ 
pression  cay  be  stated  as 


-  (( 


&  „ (  or  (Ai)))t 

k<Ri  jfR'fk  J 


’  S^) , 


(4) 


where  che  arrow  indicates  furry  sec  icplicacion.  Reducing 
che  icplicacion  to  core  pri-ative  futty  sec  operations 
and  enploying  futty  set  ceabership  function  notation, 

One  obcains 


0or(l-fo  <£or  (*Aj  f'j))).  0B.«+» 

1  k(Rl  J(llk 

-  .1  ~  •**  *  .  *  ' 

where  0  is  a  large  vector  consisting  of  all  the  Vj's  and 
is  some  attribute  delineating  the  possibility 
distribution  of  the  correlations. 


Many  of  these  isplication(or  modus  ponens)  rules, 
before  being  put  into  the  forms  given  in  eq.(4)  or  above, 
are  typically  obtained  from  the  panel  of  experts  in  an 
observed  data-decision  operational  framework,  as  the  fol¬ 
lowing  examples  shovr 


1.  If  targets  A  and  3  are  such  that  their  (observed) 
signal  characteristics  catch  reasonably  well  (this  can 

be  made  core  specific),  then  they  probably  correlate. 

2.  If  targets  A  and  B  arc  such  that  when  updated, 
their  regions  of  uncertainty  reasonably  overlap  (again, 
this  can  be  core  specifically  quantified),  then  they  arc 
candidates  for  correlating. 

3.  If  targets  A  and  B  match  on  certain  character¬ 
istics  but  not  oc  others,  then  correlating  may  or  cay 
not  occur.  (In  practise,  this  rule  would  be  replaced 
by  a  number  of  rules  containing  various  combinations  of 
catching  attributes  and  various  conclusions  as  to  che 
possible  correlation  levels.) 

4.  If  targets  A  and  3  arc  such  that  their  positions 
match  up  to  some  gating  level  C  and  their  visual  forms 
appear  to  match  up  to  some  gating  level  D  (as  for  ex¬ 
ample  by  comparing  their  lengths,  shapes,  markings, 
etc.),  then  they  most  likely  correlate. 


A  brief  comment  on  che  operators  (k  and^  is  approp¬ 
riate  here.  These  fuzzy  sec  operators  represent  'and* 
and  'or*,  respectively.  (Incidcncly,  'not*  is  repre¬ 
sented  by  the  operation  l-(-).  throughout,  but  will  not* 
be  explicitly  used.)  In  the  Initial  period  of  fuzzy  set 
theory  (circa  1965-1975),  these  operations  were  usually 
Interpreted  as  min  and  max,  respectively.  However,  more 
recently,  it  has  been  shown  from  both  an  empirical  view¬ 
point  and  theoretical  considerations  that  these  opera¬ 
tors  should  have  a  more  flexible  interpretation.  One 
class  of  such  interpretations  (theoretically  justified) 
consists  of  the  triangular  or  t-norms  and  conorms.  Thus 
the  functions  prod  and  probsum  as  well  as  many  others 
may  well  be  used  in  evaluating  these  operations.  (Sec 
Zimacrmann  (23)  for  empirical  studies  and  Klement  (24) 
and  Coodman  (11)  for  theoretical  work  in  this  area.) 
Throughout  this  paper  wc  do  not  specifically  evaluate 
the  *  and'  and  'or*  operations.  Work  is  currently  being 
carried  out  to  determine  which  evaluations  arc  cost 
appropriate  for  che  ocean  surveillance  correlation  prob¬ 
lem  and  a  future  paper  by  the  auchor  will  discuss  this. 
(Sec  RESEARCH  ISSUES.)* 

Analagous  to  the  situation  in  step  III,  each  rule  is 
actually  obtained  from  the  panel  of  experts  and  thus  cay 
vary  somewhat  from  individual  to  individual.  Thus,  in¬ 
itially,  equation  (5)  should  be  codified  to  reflect 
as  is  che  case  In  equation  (3). 

VI.  In  steps  III  and  V,  the  possibility  distributions 
for  the  attributes  and  the  rules  dcpcnd^n-esch  expert's 
intcrpretaticn.  A  procedure  is  required  which  will 
average  out  these  variations, yet  retain  all  of  the  in¬ 
formation.  Two  options  arc  available.  _Jhe  first  em¬ 
ploys  a  probabilistic  approach.  Each  expert's  response 
to  a  particular  attribute  or  rule  Is  weighted  and  the 
corresponding  possibility  distributions  are  summed 
poinevise.  The  convergence  of  this  expression  to  che 
'true*  value  can  be  justified  by  appealing  to  the  Law 
of  Large  Numbers.  (See  Coodman  (25), Theorem  2.2)  for  a 
sore  thorough  discussion  and  results.)  The  second 
option  in  combining  che  possibility  distributions  is 
to  use  a  fuzzy  set-multiple  valued  truth  theory  approach. 
In  this  approach,  each  response  by  — »  expert  is  con¬ 
sidered  a  conditional  one  and  corresponds  to  a  condi¬ 
tional  fuzzy  set  (her.ee.  che  notation  in  cq.  (3),  for 
example)  relative  to  the  Inecrprccacion  of  the  operators 

and  $or  .  By  developing  a  general  theory  for  such 

conditional  fuzzy  sets,  and,  in  turn,  deriving  an  ex¬ 
tended  fuzzy  see  Bayes* theorem,  a  fuzzy  sec  form  of  che 
Law  of  Large  Numbers  and  Central  Limit  Theorem  can  be 
ocaincd.  (Again,  see  (25) .Theorems  2.1  and  2.2)  as 
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well  as  (II),  section  3).  In  particular,  non-trivial 
results  stay  be  derived  for  the  interpretation  #or) 
■  (prod,probsu=) ,  but  Ironically,  not  for  (aln.oax)! 


C(C,g)  Is  replaced  by  any  proj(GlC.g)) ,  that  is  hny 

(fuzzy  set)  projection  into  anv  subspace  of  X- 

(Proof:  Follows  easily  fros  the  theores.  again,  sec  (25)) 


As  a  consequence  of  the  above  discussion,  froa  now 
on,  the  Individual  expert  variation  response  will  be 
oaitced. 

VII.  The  heart  of  the  procedure  utilizing  attribute 
information  is  based  on  a  deductive  logic  theorea  first 
exhibited  in  a  core  narrow  fores  in  (3),  and  then  extend¬ 
ed  to  a  very  general  setting  In  (25),  Theorea  2.3  . 

Theorea  Uniformly  cost  accurate  confidence  sets. 

Let  O  (Cj  |  j«l,..,s}  be  any  collection  of  fuzzy  sub¬ 
sets  of  sose  fixed  base  space  X.  Thus,  dC|  :  X  •*  CC,  I3t 
J*l,2,...,n  ,  are  the  corresponding  possibility  distri¬ 
bution  functions.  Let  La  ■  the  set  of  all  l  by  a  row 
vectors  a  -  (a|...,a_),  where  each  a;  is  such  that 
0<  ajil-  Let  g  be  any  function  where  g:La  -  (O.l)ls 
such  that  g  is  non-decrcaslng,  l.e.,  if  a*  ■(a^,...,aL) 

5  a*«(ai, . . .  ,aa)  chat  is,  a'i  <  a'i,  for  all  i),  then 
g(a')  $  g(«") 

Define  the  fuzzy  subset  of  X,  C(C,g)  by  the  possi¬ 
bility  distribution  function  oc^c  ,  -[0,13  ,  where 
for  any  WcX,  *5 

*C(C,g)(«)  *  g((4cl(«} . ♦ca<W»>-  (*> 

Kexc,  for  an}-  acts,  define  che  ordinary  subsec  of  X 

H(C.a)  •  |  sCj(«)  >  aj,  utx,  for  all  j  }  .  (?) 

and  for  any  fuzzy  subset  A  of  X  and  real  nuaber  u, 

0  <  u  <  1,  define  the  ordinary  subset  of  X 

-«K(a.u)  -  {-u-|  »A(u)  t  u,  «cx4  --  -  -(S).  _ 


Then: 


(i)  H(C,a)  £  K(C(C,g),g(a)) ,  for  all  C,a  as  above. 

(ii)  If  A  is  any  fuzzy  subset  of  X  such  that 
H(C,a)  S  K(A,g(a)),  for  all  a  as  above, then 
necessarily 

K(C(C.g).g(a))  S  K(A.g(a)).  for  all  C.a  as 
before.  «- 

(For  a  proof,  see  (25).)  ^ 

Thus,  K(C(C,g),g(a))  can  be  considered  to  be  che 
'tightest*  (ordinary)  g (a) -level  confidence  subset  of  X 
containing  hypotheses  set  H(C,a),  simultaneously  for  ail 
acLc  .  Exaspics  of  function  g  in  the  theorem  include: 
any  weighted  average,  prod,  sin,  max.  and  many  ocher 
functions,  including  all  t-norms.  Indeed,  it  can  be 
shown  that  if  g  is  chosen  to  be  the  c-norm  used  for  che 
interprecacion  of  the  *and*  operator  then  che  con¬ 
fidence  set  K(G(C.g) ,g(a))  also  represents  the  fuzzy 
set  interpretation  of  the  entire  hypotheses  set  H(C.a). 
On  the  other  hand,  use  of  che  weighted  average  leads  to 
core  stable  values  for  the  confidence  g(a)  than,  for 
example,  use  of  the  t-norm  min  (which  can  be  lowered 
drastically  by  merely  one  snail  confidence  value). 


Corollary 

With  the  sane  assumptions  in  che  above  theorea.  all 
results  carry  over  with  the  obvious  sodif ica cions  when 


VIII.  The  theorea  in  step  VII  is  applied  to  the  situ¬ 
ation  at  hand:  Specific  rules  are  selected  as  appropriate 
and  attribute  data  is  observed.  Thus 

(O  Zj  is  observed  with  some  confidence  o<  •  where 
any  confusable  attribute  value  Vj  for  Aj  (in  Jdomaln  Dj) 
satisfies  the  confidence  relation 


^Ajlz/V  ‘  “j  ■  W 

using  the  modeling  froa  step  III. 

(li)  Rule  Ei  is  selected,  for  say,  1-I....S.  Each 
such  rule  is  assigned  sose  confidence  level  so  chat  the 
arguments  V  and  Q+  jointly  satisfy  the  relation 


«Ei  (V.Q+)  >  3t 


(10) 


(111)  If  any  rules  or  attribute  data  is  Initially- 
described  in  terns  of  randos  confidence  sets  -  which 
are  usually  In  a  conditional  parameter  fors-a  procedure 
exists  for  first  replacing  these  confidence  sets  by 
chose  not  In  conditional  foes  and  then  converting  to 
posslbillscic  lores.  (See  (13),  Appendix  A  and  eq.(2.t)3 
(lv)  In  the  theores.  replace  each  C«  by  an  (AjIZ.)  or 
Ej: ,  whichever  are  appropriate,  and  slsllarly  replace 
each  aj  by  either  Sj  or  2^  . 

(v)  Choose  for  g  sone  convenient  function  such  as  an 
averaging  one.  sin,  or  prod.  The  appropriate  choice  for 
g  can  be  decerslncd,  for  exasple.  by  che  choice  of  c~ 
norn  for  the  ’and'  operation.  Although  weighted  averag¬ 
ing  does  not  yield  a  t-nora,  It  is  also  a  natural  way  to 
obcain  a  single  f lgure-of-scric  confidence  level  fra= 
the  given  ones  (l.e..  the  a,‘s  and  3,'s),  froc  a  statis¬ 
tical  view  point.  _(See  the  resarks  'following  the 
theores  in  step  VII.) 

(vi)  Use  the  corollary  following  the  theores  In  step 
VII.  where  the  projection  operation  Is  into  che  space  2 
of  all  possible  correlations  (or  data  partitionings} 

Q+  ,  el isinating  finally  the  ’nuisance’  vector  dosatn 
corrcspond-to  all  possible  values  for  the  0. 


Specifically,  subsceps  (i)-(v)  above  lead  to  the 
fuzzy  sec  G(C,g)  which  Is  described  by^csslblllcv  dis¬ 
tribution  s)  described  In  eq.  (6 Tr,  where  argusent 

»-(V.Q+) .  Then  proj  (C(C,g))  here  is  decersined  by  the 
relation 


Voj-(C(C.g))(Q+>  "  *or  <sC(C,g)(V’Q+,)’  Cl  I) 
*  Vc  V 

where 2/ la  jhe  sec  of  all  possible  v's,  I.C..V  ■ 

u 

^  (Dj),  and  $  is  sose  appropriately  chosen  t-conota 
j-I 

for  che  fuzzy  sec  systes.  (See  the  consents  near  che 
end  of  step  v.) 


The  final  result  of  the  above  eight  steps  is  (via 
che  corollary  of  step  VII)  that  Kfproj  .^(CtC.gD.gfa)) 
Is  che  tightest  or  uniforsly  =o«c  accurate  g(a)-le“el 
confidence  subset  of  2.  eontaining  the  hypotheses  sec 
projection  into  2  ,  proj  (H(C.a>).  Thus,  in  an  opcisal 
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nathematical  way,  this  confidence  sec  describes  the 
possible  values  for  given  all  of  the  relevant  at¬ 
tribute  information  and  rules.  (The  following  example, 
gives  a  sinple  geometric  illustration  of  che  above.) 


UNCLASSIFIED 


Conclusion  1C.  trftfi  boundary  contains  at  a  propar  suD**t.  prerrj (At  A  A2). 
Indicated  by  '/rf/.  Projection  *C,»  PROJ*  (KJ  it  Indicated  by ...... 


(U)Ccoaccrlc  Interpretation  of  Conjunctive  Premises 
Figure  5 


The  second  procedure  requires  additional  imputa¬ 
tions.  but  is  based  on  core  rigorous  grounds. 

Recall  that  is  the  class  of  all  possible  correla¬ 
tions  Q+  .  Let  ^(2)  denote  the  class  of  all  fusty 
subsets  A  of  3  with  corresponding  possibility  distri¬ 
bution  function  $A  :  3  -  CO.ll  .  Note  chat  Aq  s 
pr°ja(C(C.g})  cA{3)«Lec  A(2  )  denote  che  class  of  all 
random  subsets  of  2  .  (See  (13)  for  background.)  Then 
it  is  knoun  ((10). (11))  that  a  number  of  such  cappings 
S  exist,  called  choice  functions,  such  that 

S:  /=d)  ~/?(2)  (onto)  .  (13) 

where,  for  any  A  t  ^O).  S(A)  t  /?(2)  Is  such  that, 
for  all  Q+  c  2  , 


P<Q+  c  S(a))  -  4A(Q+) _  (!l 

and  such  chat  certain  homomorphic  fuzzy  set  and  random 
set  operator  relations  hold  for  S  and  a  related  fa=ily 
of  cuppings.  (See  (10)  and  (11).)  Each  S  generates  the 
equivalence  class  of  random  sets  mentioned  briefly  in 
the  Introduction. 

Two  of  the  cost  important  choice  functions  relating 
fuzzy  sec  theory  and  probability  theory  as  outlined 
above  are: 

(l)  S  ■  Su  .  where  0  is  a  randos  variable  distri¬ 
buted  uniformly  over  CO.IJ  .  and  for  any  A  c  ^(S), 


Clearly,  this  procedure  can  be  expanded  directly  to 
Include  all  Information  concerning  correlation:  in  parti¬ 
cular,  geolocational  information  can  be  easily  Integrated 
Into  the  scheme.  However,  as  stated  at  the  outset  of 
this  paper,  since  geolocational  and  false  alarm  inform¬ 
ation  has  been  treated  historically  from  a  different 
viewpoint-.- namely,  a- probabilistic  approach-,  which  al«o  - 
has  many  justifications  for  us.  ,  an  all-encompassing 
fuzzy  set  approach  to  che  contact  c-  relation  will  not  be 
persued  here.  Future  work,  however,  may  lead  at  least 
partially  In  this  direction. 


%l(A)  -  »_1  (CU.il)  -  {Q.  |  4A(Q.)  i  u}  (15) 

and 

(il)  S  •  T.  where,  for  any  A  tA3),  T(A)  is  de¬ 
termined  by  considering  the  corresponding  (ordinary 
random)  membership  function  where  each  of  the 

latter's  one  dimensional  zero-one  marginal  random  vari¬ 
ables  *j(a)«J.).  Q+  t  2  ,  are  all  mutually  statistically 
Independent  with 


DETERMINATION  OF  POSTERIOR  POSSIBILITY  ASD  PROBABILITY 

'distributions 

Ulch  the  determination  of  the  optical  confidence 
sec  for  Qt  ,  given  all  attribute  data  and  rules,  the 
posterior  possibility  distribution  function  for  0^  can 
Chen  be  determined.  Clearly,  by  Inspection  of  the  form 
of  che  optimal  confidence  sec,  Ic  follows  that  the  de¬ 
sired  distribution  function  Is  given  by  the  truncated 
form 


4?(Q„)  .  |VoJa<e{C.g)>«4>.  H  0*  t  Ko 

l  0  .  if  Q*  ♦  So  , 

where  K„  a  K(proja(C(C.g)),g(a)),  noting  that 
0+  t  So  iff  *proja(G(C,g))(<l+)  *  s(,)- 


(12) 


!P(«T(A)(Q+>  -  1)  -  p(Q+.  C  T(A))  -  J .«)*.> 

(16) 

S><«T(a)M*>  *  0)  -  p(Qa  t  T(A))  *J-»A«+>  . 

J  * 

Suppose  now  chat,  -ithout  loss  or^geaerallty. 

*o  •  . .  "here 

0<*U)£4Ao<^(1>)spA3(<ll.<2')<...<pAoCQt{::,>^«l  .  (I») 


(The  following  results  can  be  modified  if  equalicy 
occurs  in  places  In  eq.  (17).)  If  S  is  any  choice 
function,  there  is  a  naturally  corresponding  random 
variable  IX  (SfA,,))  over  3.  representing  the  maximal 
possible  element  of  random  sec  S(A„)  with  respect  to 
So  given  by 


Decerminacion  of  che  posterior  probability  distri¬ 
bution  function  corresponding  to  4p  can  be  accomplished 
In  two  different  ways. 

The  first  procedure  is  an  immediate  consequence  of 
eq.  (12):  Either  Identify  directly  or  renormalize  che 
possibility  distribution  ip. 


K(S)Aq)}  -*  that  unique  value  of  Q,,  i.e«,  some 

for  ISjim  such  that  maxfc^Q^))  , 

f(S(Ao)nEo)J 

occurs  for  <L<J>  •  provided  that 
S(Ao)fl  Sy,  7  <ZS  .  (IS) 
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and 

It  (S(A<,))  "  5*  “hen  S(*o)  f)  K„  -  &  ,  (19) 

Ic  folio VS  chat  the  probability  function  for 
is  given  by,  for  any  j,  Ijjia  . 

?at(s(Ao)) .  •  "£  p<s(a0)  .a)  . 

a  cj-  (20) 

S-  -  {C?  |  <X  CK,,  and  Z?(C?)  ■  . 

Thus,  the  following  idcncif ication  csay  be  cade: 

pW+I  2a)  *  P(«(S(A0))  -  Q+)  .  (21) 

In  particular,  specialization  o £  the  above  results 
to  cases  (1)  and  (11)  for  choice  function  S  yields: 

(111)  For  S  -  S0  , 


pfWfSfAo))  -  <),.<=>)  -  4A  (Q>>)  . 

o 

(22) 

p(7f(S(A0))  -  Q+O))  .  0  .  lijSa-1  . 

(23) 

p(Z«S(Ao)>  -  f#  )  -  l-:j^Q,(=>>  . 

(24) 

i  For  S  -  T  . 

p(2?(S(Ao))  -  Q+U>) 

*A0M+<i>)  *  7jlO-4Ao«4(i>))  • 

(25) 

for  Osjsm  ,  where  if  j  ■  n.  the  product  term  In  eq. (25) 
Is  defined  to  be  unity,  and  if  j  -  0.  Q,.*0)  Is  defined 
to  be  equal  to  the  null  si  0  . 

-  .i.  -  .  —  —  -  ,  _ 

Clearly,  use  of  T  leads  to  a  core  tractable  result- 
than  the  one  point  nass  result  due  to  Sy-  However,  for 
either  choice  function,  the  resulting  evaluation  for 
saxinal  p(Qt  I  Za)  coincide,  the  desired  value  occurring 
at  Q+  -  Q+(=),  the  cost  possible  value  of  Q^.  . 

Finally,  cabulate  Che  possibility  distribution 
Cy(Q-)  versus  Q+  (  by  first  tabulating  dA(Q*)  versus  Q+) 
for  all  feasible  Q*.  t  “X  ,  and  substitute  0  these  values 
Into  eqs.(22)-(25lC  and  using  eq.(Zl).  finally  Into 
eq.(l). 

RESEARCH  ISSUES 

Throughout  this  paper.  It  has  been  enphaslted  that 
for  any  reasonable  determination  of  operator  pair 
(Og,vor)  ,  the  entire  procedure  renains  valid,  con¬ 
ditional  upon  the  interpretation  of  these  operators, 
such  as  (sin, sax)  or  (prod.probsus).  Yet  in  order  to  is- 
plesent  the  sehese,  a  specific  evaluation  Is  obviously 
needed  for  the  'and*  and  'or'  operators- 

It  has  been  shown  in  (11), and  mentioned  briefly  In 
step  ’  of  the  basic  codeling, chat  a  reasonable  family  of 
furry  set  operators  (justified  by,  e.g. .  relations  with 
multiple-valued  logic  and  set  theory- see  (24) )  to  con¬ 
sider  Is  that  of  pairs  of  t-norms  and  c-eonoms  (usually 
restricted  to  have  DeMorgan’s  property-  see  (Il».  (See" 
(12)  and  (24)  for  background  and  further  properties.) 
Furthermore,  it  has  been  shown  in  (il)  that  a  particular 
subfamily  of  such  operators  yields  especially  close 
(homomorphic)  relations  between  all  furry  sec  systems 


determined  by  these  operations  and  corresponding  random 
sec  systems.  This  subfamily  cf  pairs  of  operations  (for 
'and'  and  *or‘)  consists  essentially  of  (min,  max), (prod, 
probsum),  and  all  countably  Infinite  or  finite  convex 
weighted  sums  of  the  above  operations  restricted  to  dis¬ 
joint  regions  (chough  somewhat  restricted  In  form), 
which  include  the  first  two  pairs  as  special  cases.  Cur¬ 
rently.  work  Is  going  on  In  determining  empirically  what 
values  the  above-mentioned  weights  should  be  assigned. 
(This  Is  somewhat  analagous  to  the  work  of  Zimmermann 
(2)),  who  used  a  different  family  of  operators  and  empir¬ 
ically  estimated  certain  adjustable  parameters  In  the 
family.) 

Alternatively,  another  family  of  operators,  pos¬ 
sessing  some  (but  not  all)  of  the  desirable  properties 
of  the  first-mentioned  family,  has  been  shown  to  yield 
very  desirable  furry  set  analogues  of  the  Laws  of  Large 
Numbers  and  Central  Limit  Theorem.  (This  family  Is 
denoted  as  the  Archlmedlan  Franklan  family  and  is  dis¬ 
cussed  In  (11).)  On  the  other  hand,  (sin,  max)  does  not 
yield  desirable  asymptotic  system  properties  (although 
(prod,  probsum)  does,  as  an  internal  member  of  the  Arehi- 
medlaa  Franklan  family),  r.or  apparently  do  any  of  the 
weighted  sjms  previously  referred  to.  A  paper  on  this 
topic  will  be  forthcoming. 

Trade-offs  need  be  established  between  utilities  of 
choice  of  the  various  plausible  operator  pairs  for 
^1  •  Air)  •  before  determining  a  final  candidate  pair. 

SVXMASY 

Am  approach  to  the  utilisation  of  attribute  Informa¬ 
tion  for  che  contact  correlation  problem  has  been  out¬ 
lined  In  this  paper.  The  novelty  of  the  technique  lies 
in  the  use  of  possibility  theory  In  the  modeling. 

The  required  steps  in  developing  this  technique 

were: 

(I)  Establishment  of  a  set  of  relatively  primitive 
attributes. 

(II)  Determination  of  attribute  domains 

(III)  Querying  of  a  panel  of  experts  to  establish  post¬ 
erior  possibility  distributions  directly.  In 
place  of  a  Bayesian  approach,  for  the  accrlbute 
values. 

(IV)  Reducing  che  calculation  load-i.y_use  of  cercain 
analytic  models  such  as  thosa^if  the  translation 
type. 

(v)  Determination  of  modus  ponens  rules  which  delin¬ 
eate  the  possible  correlations',  by  again  ques¬ 
tioning  the  available  experienced  personnel. 

(V!)  Smoothing  out  of  the  variability  in  the  models  ob¬ 
tained  In  steps  (III)  and  (V)  due  to  individual 
responses.  Either  probabilistic  or  analagous 
furry  set  asymptotic  results  are  deployed. 

(VII)  Demonstration  of  a  general  theorem  which  shows 
that  given  a  set  of  hypotheses  H  formed  by  the 
conjunction  of  individual  confidence  sets  (des¬ 
cribed  by  possibility,  or,  in  effect,  probability 
distributions),  a  uniformly  most  accurate  confi¬ 
dence  set  S  exists  described  by  a  single  possi¬ 
bility  distribution,  which  contains  H.  The  latter 
possibility  distribution  was  shown  to  be  con¬ 
structed  by  a  simple  application  of  any  one  of  a 
large  class  of  (non-decreasing)  functions  which 
also  determine  the  tonfidence  level  for  I. 
Modification  of  this  theorem  for  projection  oper¬ 
ations  was  displayed  in  the  form  of  a  corollary. 
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(VIII)  Application  of  Che  corollary  given  In  seep  (VII) 
co  the  attribute  problem.  assumed  to  be  in  the 
foes  of  a  conjunction  of  attribute  data  and 
available  nodus  ponens  rules  (nodeled  In  step 
(V)). 

(IX)  Derivation  ot  posterior  possibility  distribution 
for  the  correlations  given  the  attribute  informa¬ 
tion.  In  turn,  the  posterior  probability  distri¬ 
bution  function  for  the  correlations  vas  also 
obtained  by  two  different  approaches,  the  first 
being  relatively  slnple,  the  second,  requiring 
sore  computations,  but  derivable  from  a  sounder 
basis  (using  relations  between  fuzzy  and  random 
sees). 

(X)  Substitution  of  the  results  from  step  (IX)  into 
the  basic  factor  model  for  the  overall  posterior 
probability  distribution  of  the  correlations. 

Implementation  of  the  technique  demands  a  specific 
choice  for  the  'and'  and  'or'  fuzzy  set  operations  used 
throughout.  Some  discussion  vas  presented  concerning 
this  problem. 

The  same  procedure  developed  in  this  paper  could 
also  be  used  in  a  vide  variety  of  problems  involving  the 
estimation  of  an  unknown  parameter  vhen  some  of  the 
available  information  is  given  in  linguistic  form.  (See 
some  earlier  related  uork  in  (3),  especially  section  3.) 

Future  work  will  be  greatly  concerned  with  real- 
vorld  implementations  and  modifications  of  the  technique 
presented  here.  The  repotted  successful  implementation 
of  a  number  of  fuzzy  set-logical  approaches  to  non-mili¬ 
tary  problem  areas  such  as  in  medical  diagnosis  (25)," 
library  search  syscens  (37).  and  fault  analysis  (33). 
nay  veil  serve  as  an  impetus  fot  the  treatment  of 
military  problems  by  such  techniques,  as  presented  in 
this  paper. 

.  L  -  — .  “  -  . 
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